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'~{ I We elucidate the universal properties of the nonequilibrium steady states (NESS) in a driven 

symmetric binary fluid mixture, an example of active advection, in its miscible phase. We use the 
symmetries of the equations of motion to establish the appropriate form of the structure functions 
OA ' which characterise the statistical properties of the NESS of a driven symmetric binary fluid mixture. 

OA , We elucidate the universal properties described by the scaling exponents and the amplitude ratios. 

Our results suggest that these exponents and amplitude ratios vary continuously with the degree 
of crosscorrelations between the velocity and the gradient of the concentration fields. Furthermore, 
we demonstrate, in agreement with Celani et al, [Phys. Rev. Lett., 89, 234502 (2002)], that the 
^ . conventional structure functions as used in passive scalar turbulence studies exhibit only simple 

scaling in the problem of symmetric binary fluid mixture even in the weak concentration limit. We 
also discuss possible experimental veriflcations of our results. 
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, Critical dynamics of equilibrium systems show universal scaling properties for thermodynamic functions and corre- 
lation functions characterised by universal scaling exponents that depend on the spatial dimension d and the symmetry 
I ' of the order parameter (e.g., Ising, XY etc.) and conservation laws 0,01) but not on the parameters that specify the 
' O details of the Hamiltonian. A different situation arises in driven, dissipative, nonequilibrium systems with noncquilib- 
C ' rium statistical steady states (NESS). In this letter we demonstrate that in the NESS of externally driven turbulent 
homogeneous and isotropic symmetric binary fluid mixtures (50-50 concentration of the two components) universal 
properties are characterised by certain structure functions, determined by the symmetries of the equations of motion, 
and which arc different from those used in pasive scalar turbulence studies. We find that the universal properties vary 
'T'* ' continuously with the degree of crosscorrelations between the velocity and the gradient of the concentration fields, a 
^ [ situation which has no analogue in equilibrium. 

^y-^ ' The turbulent hydrodynamics of an incompressible symmetric binary fluid is described by the coupled evolution of 
ly-j ^ the velocity field u(x, t) and the concentration field gradient VV'(x, t) with ^/'(x, t) = t) — pb{^, 0]/Po where po 

. is the mean density, and pA and pB are the densities of the two components. In a symmetric binary fluid t)) — 0. 
' Experimentsperformed on such systems have been typically concerned with the measurements of effective transport 
coefficients |3|. In a symmetric binary mixture, ip is not simply advected passively by the u- field, but is active, i.e., 
■ concentration gradient Vi/j reacts back on u to produce flows. In ordinary fluid turbulence the NESS is characterised 
by the multiscaling exponents of the structure functions of the velocity fields (see below) H Compared to the 
studies of ordinary fluid turbulence the study of turbulent binary mixtures is in its infancy [Sj . 

In this letter we develop a one-loop, self-consistent theory for a randomly forced symmetric binary mixture system 
. to elucidate the universal statistical properties of its NESS in the miscible region. We identify the symmetries of the 
(~| Eqs. of motion which allow us to construct the structure functions (different from the ones discussed in the literature 
O : for passive scalar turbulence studies - see below) which exhibits multiscaling in the NESS in the inertial range.. 
O ' We calculate the interrelations between certain dimensionless numbers, e.g., the amplitude-ratio of the equal time 
correlators of u and Wip fields and the renormalized Prandtl number (see below) in the NESS in d-dimensions. Finally, 
we elucidate the nature of multiscaling by calculating the multiscaling exponents for the even order structure functions 
of the VV' field in the limit of small ip in two different models. In one case wc find the multiscaling exponents to vary 
I with the degree of crosscorrelations between the u and VV' fields. We also consider the scaling of the conventional 
■ " ' structure functions and show, in agreement with Ref.|^, that they exhibit only simple scaling. 

a. Equations of motion, their symmetries and structure functions:- It is well-known 0il3l3 that in the dynamics 
of a symmetric binary fiuid mixture the velocity fields u couple with the concentration gradient Vip and not with 
itself. Thus it is useful to write the coupled evolution equations in terms of u and V?/'. We are interested in the 
isopropic and homogeneous situation, i.e., we have no mean concentration gradient, (V?/') = 0. Furthermore, for 
calculational convenience (see below), we write = b(x, i) -|-Bo/q!o, where b is a space-time dependent vector, Bq 
is a constant vector and ag is a numerical constant. Note that since we impose (VV') = we have (b) — — Bo/ao- 
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Such a parametrization of Wip yields the Eqs. of motion 

+ Ai(u • V)u = - Bo • Vb - aobV • b + vV'^n + f , (1) 

at po 

c^b A 

and — + — V(u • Bo) + A2V(u • b) = r/V^b + g. (2) 

ot ao 

We also consider, for simplicity, the velocity field u to be incompressible enforced by the condition V • u = 0. The 
parameters Ai, A2 and ao are coupling constants. Coupling constants Ai, A2 are unity and kept for formal bookkeeping 
purposes. A choice of ao depends on the unit of Vip and is given by ao — ^ where C is the thermal correlation 
length, X is the susceptibility per particle, and m is the particle mass. This, within the mean field theory, yields ao 
to be approximately temperature independent near the critical point Tc 7]. On the other hand, away from Tc, similar 
considerations yield ao ~ T. Therefore, at any finite temperature ao ~ 0{1) and, as a result, a symmetric binary 
mixture exhibits active behaviour. Functions f and g are stochastic forces required to maintain a statistical steady 
state. Parameters I'o, rjo are the bare viscosity and diffusivity respectively. Note that in a problem of true passive 
advection ao is identically zero. At this point we would like to clarify that by writing \/ip as the sum of b + Bo/cto 
we do not introduce any anisotropy in the problem as the latter is decided by (Vip) = (b) + Bo/ao which is kept zero 
in our system. Furthermore, stochastic forces f and g have isotropic variances [see Eqs.©] and hence do not violate 
the isotropicity of the system. 

The Eqs. of motion ^ and Q are invariant under 

• TI: The Galilean transformation u(x, t) u(x + Uo^, i) + Uo, ^ — > ^ — Uo • V and b ^ b with Ai = A2 = 1 

Hi 

• TII: The transformation Bo Bo + aoS, b(x, t) b(x, i) — 5, u — > u. Here the shift 5 is a vector (see, e.g., 
Ref.Q for an analogous symmetry in a related problem). 

Transformation TII implies that the rcnormalised versions of the Eqs. and ^ must be expressible, as the bare 
equations themselves, in terms of the total gradient fields (= b + Bo/ao) |9|. It has been overlooked in related 
previous works The existence of TI implies nonrenormalisation of Ai, A2. Since the renormalized versions of 

Eqs. (Q) and Q are are also invariant under the transformatio TII, the fourth and the fifth terms in Eq. ^ must 
renormalizc in the same way. Similarly, the second and the third terms in Eq. lI2Jl must renormalizc in the same way. 
Hence the corresponding renormalization-Z factors must be equal leading to Zb^Zi, — Z^^Z^. It can be argued that 
in model Eqs. iQJ and ^ Zb^Zi, = 1 ^3 (this essentially is a freedom of choice for the renormalization Z- factors 
in the present model). This yields Zc^^ = 1. Here, Z^^, Zb„ and Zfy are renormalization factors for ao, Bo and b 
respectively. Therefore, ao does not renormalize and hence it can be set to unity by treating all the gradient fields as 
effective gradient fields ^/a^h. In fact, this can be argued in a simple way: Let us assume that under mode elimination 
and rescaling ao Pao where /3, a scale factor, can be absorbed by redefining the units of the gradient fields b by 
^/]3h b. This leaves Eq. (j^J and ao in Eq. |^ unchanged. Furthermore, the invariance of the Eqs. motion 
and ^ under TII ensures that the composite operators (&i)" (Cartesian indices i contracted or free) are also not 
renormalized . It should be noted that these results hold regardless of the strength of the feedback of b on 

u. Therefore, even though the dynamics of the b fields are linear in its passive limit, i.e., when the feedback of b 
on u is ignored, it must still be consistent with TII and its consequences 10] for it to represent the physical weak 
scalar limit of the binary mixture problem. We emphasize that by the transformation TII we do not introduce a mean 
concentration gradient {yip) in the system as the latter, being given by Bo + (b(x, i)) is independent of such shifts. 
In particular we work with (VV') =0 (i.e., no mean concentration gradient). This, furthermore, can be achieved in 
a perturbative calculation by (i) (b(x, t)) = —Bo 7^ 0, or, (ii) (b(x, t)) = —Bo — 0, both of which are in agreement 
with (VV') — 0. The invariance under the transformation TII ensures that physical results are independent of such 
choices. In our calculations below we work with (b(x, t)) = —Bo = which is the most convenient one. 

In fluid turbulence, the NESS is characterized by the multiscafing of the structure functions S^{r) — {[\ui{x + r) — 
Mi(x)]fi]") ~ r^" for r in the inertial range (i.e., rjo ^ t <^ rjo is the dissipative cutoff scale, L is the integral scale 
~ system size) which are invariant under the Galilean transformation. Analogously, in the problem of symmetric 
binary fluid mixtures, the appropriate structure functions must be invariant under TI and TII. Such functions are 
S'°(r) = ([|ai(x + r) — ai(x)]fi]"), a = u,b. The scalar field structure functions as used in passive scalar turbulence 
studies [Tl|, 5*^ = {[ip{x + r) — '(/;(x)]") are not invariant under TII: b ^ b + constant, and exhibit only simple 
scaling in the symmetric binary mixture problem ( and hence not suitable for the studies of the passive limit of 
the active problem), as we show below and has been observed in Ref. 0. In this communication we calculate the 
multiscaling exponents of the even order structure functions S2n{r) = ([&i(x + r) — 6i(x)]^") by using one- loop dynamic 
renormalisation group (RG) in a field theoretic framework in the weak scalar limit of the symmetric binary mixture 
problem. 
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b. Correlation functions, scaling exponents and dimensionless numbers:- Since u and b are polar and axial vectors 
(Note: a mass density is a pseudo scalar in Sd and hence b — Vip is an axial vector) respectively, the crosscorrelation 
tensor Cf^ = (wi(k, k, 0)) is imaginary and odd in the wavevector k. The propagators Ga, a = u, 6 of u and b 

as functions of k and frequency w are defined as 

Ui(k, uj) = Gu{k, uj)fi{k, uj), bi{k, uj) = Gb{k, w)gi(k, w), (3) 

and the correlators 

Ct^{k,u) = (w,(k,c.)u,(-k,-u;)), C\^{k,u) = (6,;(k,c.)6,(-k,-^)). (4) 

In the scaling limit, in terms of the dynamic exponent z and the roughness exponents x„, Xh of the fields u and 
b, G^'^ — —iu) — G^^ — —iuj — Sb(fc,w) where self-energies Su = k~^mu{uj/k^), Sb = k^^mb{cj/k^) and 

the correlators C^- = P,j{k)k-'^-^^^~^mi{Lj/k''), C^^ = Q,y (k)fc-'^-2xb-^TO2(w/fc^) with Py(k) = S,j - ^ and 
= -pi being the transverse and longitudinal projection operators respectively. The Ward identities arising from 
TI and Til yield Xu = Xb = X ^-nd X + ^ = 1- 

We assume zero-mean Gaussian distribution for f and g with variances 

(/i(k,i)/,(-k,0)) = 2D{^{k)d{t), 
(5i(k,t)5,(-k,0)) = 2Df^{k)S{t), 

{Mk,t)g,i-k,0)) ^ 2^Df\k)Sit), (5) 



where D(j^{k) = Hij{k)k^ ^ e > 0, Hij = Pij,Qij for f,g respectively. Specifically, e — A reproduces 
Kolmogorov-type (hereafter K41) scaling for the kinetic energy spectrum [k'^~^ {u{k, t) ■ u(— k, t))] and the gradient 
energy spectrum [A:'^^^(b(k, i) • b(— k, i))] Allowing 1?^^ (k) , (odd in k), to break rotational invariance or 

symmetry with respect to an interchange of the cartesian indices i and j we write D^j (k) = D^j ^ (k) + D^j ^ (k) , 
with Df-^^^ and Df-'^^'^ being the symmetric and the antisymmetric parts of Df-''^\ We choose Df-'^'^'' (k) = 

D°Pij{)<i)k^~'^~'^, Dij''^^ (k) — D°eijpkpk~'^ , eijp is the totally antisymmetric tensor in 3c? 15]. A one-loop self- 
consistent mode coupling (SCMC) theory is conveniently written in terms of the self-energies and the correlators 
involved. 

We start with the zero frequency {uj — 0) forms for the self-energies 

E„(fc,tj = 0) = i/F, Sfc(A:,w = 0) = r7fc^ (6) 



and the correlators 



2 D 

G^^{k,^ = Q) = --^p^k-^^-d-\ 



77- 
1 

vrj 



Gt^{k,uj = 0) = '^^I^P,^k-^x-d-z^ JJ^^y^ ^ -Ds{~k) (7) 



in d dimensions and 



^ij — j^j^ ^ijpkpk 2X 2 4 ^g-j 

in 3d. Subscripts s and a refer to the symmetric and antisymmetric parts of C^^ . The lack of vertex renormalizations 
in the present problem in the zero wavevector limit implies that the SCMC yields exact relations between x and z, 
viz, X + 2 = 1 (see, e.g., By matching the self-energies and the correlation functions calculated at the one-loop 

order at a; = with their zero frequency expressions above (we take Ai = A2 = = 1 without any loss of generality) 
we obtain 

z = 2-e/3,x = -l + e/3. (9) 
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Furthermore, we obtain self-consistent relations between the two dimensionless numbers, the amplitude-ratio T = 
D2/D1 and the renormalized Prandtl number P = 77/^. We consider the cases of symmetric and antisymmetric 
crosscorrelations separately When Da = 0, i.e., in the absence of antisymmetric crosscorrelations, we find 



d- 


1 


1 r 


1 2(<i- 




^ 2P^((i+2) 


d- 


-1 


r(d-i) 



(10) 



d{l + P) Pd{l+P) 



d^-2 I _ 4(rf"-2)jV; 

p _ d(d+2) P'^d(d+2) P(l+P) 

^ 2r-i(rf-l) I 2jVf(d-l) ' ^ > 

p{i+p)d + P(i+P)r2d 



where 



(M,(k,t)fc.(-k,t)) "' 2 
(ui(k,t)Mi(-k,t)) 



parametrising the strength of symmetric crosscorrelations in the system. 

Similarly, for finite Da, i-G-, finite antisymmetric crosscorrelations and with Dg = 0, Eq. (|10|) remains unchanged 
(meaningful in 3d only) but Eq. changes to 



15 iSW 3P(1+P) 
4r-i 4jVg 
3P(1 + P) ~ P(l+P) 



(12) 



with A'^a = 



(Mi(k,t)bj(~k,t))„(Mi(k,t)bj(-k,t))„ 
{«i(k,t)«i(-k,t))^ 



parametrising the strength of antisymmetric crosscorrelations in the 

system (subscript a refers to the antisymmetric part). Note that both and Na are dimensionless numbers since 
Xu — Xb- Our results imply that in the NESS, the values of V and P (i) are related to each other and (ii)they depend 
on the dimensionless numbers Ng^ Na, i.e., the degree of the crosscorrelations between u and b. It is significant 
that similar variations of the dimensionless numbers and their relations with the breakdown of the invariances under 
rotation and the interchange of the cartesian indices, are observed also in MHD turbulence This strongly suggests 
that such properties, having no equilibrium analogue, are rather generic properties of a class of nonequilibrium systems. 
We close this Section by noting that the parameters Ng and Na being the dimensionless ratios of aprropriate correlators 
in the NESS are measurable in experiments and numerical simulations, and not the same as the amplitudes of noise 
crosscorrelations (jSJ which are introduced in our framework as convenient means to calculate different correlation 
functions. 

c. Multiscaling exponents C,2n'~ ^ ^^^1 characterisation of the NESS for a turbulent binary mixture requires the 
calculations of <S'°(r), a = u,b (see above) which scale as ~ r'»" for r in the inertial range. There is no controlled 
perturbative calculation which permits to do this to this date. Complications of the full nonlinear problem initiated 
studies of multiscaling for the gradient field in the passive advection limit where the backreactions of the gradient 
field on the dynamics of the velocity field are neglected. There is a large body of literature on this problem, known 
as the passive scalar turbulence, see, e.g., [lH. However we differ from others in ensuring that our calculations obey 
the symmetry TII which was not considered in the existing works. 

K41-type arguments for fluid turbulence |l3j | can be easily extended to a binary mixture system by demanding that 
the scaling in the inertial range is determined by the local length scale r and the rate of mean dissipation of energy 
(kinetic or gradient). This yields, = n/3 as in fluid turbulence. However, departures from the KAl scaling (or 
multiscaling) in fluid turbulence suggest the presence of similar features in the binary mixture problem. Below we 
analyse the scaling properties of the even order structure functions 52„(r) — ([6i(x -I- r) — 6,i(x)]2") in the passive limit 
of the full nonlinear problem. Further, the velocity field u is assumed to obey a Gaussian distribution with zero mean 
and a variance llj 

(u.(k,i).,(-k,0)> = jj^,f^^l§^6{t), (13) 

where < e < 2. The field b is then governed by the advective equation 

Bh 

— + A2V(u.b) =77oV2b + T. (14) 
ot 

Here T is a short range force with correlations (^^(x -I- r, t)Tj{x)) = QijA{mr)S{t) with m ~ 1/L, L being the system 
size and A{r/L) const, for L — > 00 Note that following the formulations of Ref.lll[ the longrange stochastic 
forcing function g in Eq. Q is replaced by the short ranged force T in Eq. (fTH) . We continue to impose the 
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consequences arising from TIL Thus it can be called the modified Kraichnan model, modified to study the weak 
concentration gradient limit of the active problem. 

We employ a one-loop dynamic renormalisation group analysis with a minimal subtraction scheme in our treatment 
[T^. Following standard procedures [Tol we obtain the critical dimension of the fields b to be given by 

where e is defined through the relation H13() . Hence, from the definition of critical dimensions 
we find that the structure functions iS2„ scale as 52„(r) = r"'^^~^^/„(mr), with yet undetermined scaling functions 
fniiTir). Anomalous scaling or multiscaling would arise if the scaling functions were to behave in the asymptotic limit 
r/L ^ as [r/LY'^ (for (5„ — simple scaling would prevail). Ignoring such possibilities for the time being we find 
that C2„(if41) = n{2 — e) which is the analog of K41 scaling in fluid turbulence in this problem [loj |. 

We now turn to the study of the asymptotic forms of the scaling functions for mr 1, i.e., in the inertial 
range rjD r <^ L by using renormalisation of the most relevant composite operators in Wilson Operator Product 
Expansions (OPE) 16] of the scaling functions. This method has already been used in critical phenomena |0, 
fluid turbulence J. 4] and passive scalar turbulence The anomalous dimensions 7„ of the most relevant compose 
operators for 0„ = [fei(x + r) — fei(x)]^" are calculated in a one-loop e-expansion. This then provides the deviation from 
simple scaling: 52„(r) ~ ^n(2-e)-i-T„ -^j^gj-g j-^ jg ^jjg anomalous dimension of 0„ coming from the composite operators. 
In the OPE of 0„ only renormalised operators which are invariant under TI and TII (which keep 0„ invariant) can 
appear. Moreover, in the passive limit the linearity of Eq. (0) ensures that the number of b fields appearing in the 
OPE of On cannot exceed 2n For O2 the leading operator is [dibj{-x)dibj{'x)] with an anomalous dimension 72 = e, 
leading to C| = 2 > C^{K4,l). For n > 1 the most relevant operator in the OPE of [9^6^ (x)]p[(9j6j (x')]*, p + q < 2n, 
which determines the anomalous dimension of 0„, is (dibjdibj)"- . This yields 0| 

CL = n(2 -a)- < C^K^l). (15) 



Note that our results satisfy key characteristics of multiscaling as observed in fluid turbulence studies |j, Q : C,2n 
a convex function of n and C2 > CtC-^^l), C,2n < C2n(^41), n > 1. The structure functions for the velocity fields S*" 
do not multiscale in the present approximation due to the assumed Gaussian distribution of the velocity fields. 

In the above we assumed the velocity fields to be Gaussian distributed. In contrast, real velocity fields follow Eq. ^ 
which we now consider. We still ignore the backreaction of b on the dynamics of u; nevertheless, our calculations are 
consistent with the requirements imposed by TII. We are able to obtain additional results concerning the multiscaling 
properties of the structure functions despite the limitations [Tol IT^ of a one-loop e-expansion method in this problem 
where e is now defined through the noise correlations (O. We follow the calculational scheme outlined above and 
elucidate the multiscaling properties of iSg^ in a one-loop e-expansion. In the present problem the critical dimension 
of the field b is (2 — e/3)/2, yielding for the structure functions 

- r"(2-/3)/„(^^). (16) 

Here e = 4 is the physically interesting case. In order to find out any deviation from this simple scaling we examine 
the asymptotic form of fnimr) in the procedure outlined above which yields the following results; Writing Sl^ni^') ^ 
^»(2-£/3)+r„ obtain pii] 

e . 4P2 ^ 

(17) 



3P{1 + P)' H-P2J 
and for n > 1, 

n{n-l)e d-1 U0P^{4: ~ d){d - 1) 



< (18) 



3P{1 + P)dU + 2 {d + 2){d + A){d + 6) 

making (^2n < C2n(^41), where P is the renormalised Prandtl number. Thus, again, exponents ^2n follow the same 
key characteristics of multiscaling as observed in fluid turbulence. The most striking aspect of our results for C|„ is 
their dependence on P. According to our results obtained above P is not a fixed number in the NESS of a binary 
fluid mixture, but rather parametrised by the degree of crosscorrelations between u and b. Even though our results 
for C,2n ^-re obtained in the passive limit their dependences on P are likely to remain in the full nonlinear problem, 
and thereby open up the intriguing possibility of variable multiscaling parametrised by the degree of crosscorrelations 
between u and b. 

Finally we consider the scaling of the structure functions ~ (['/'(x + r) — ^'(x)]^"), commonly discussed in the 
passive scalar literature We have S'^„(r) ^ r"*^'*~'^V„(TOr) in the case where the u fields are assumed to have 
a Gaussian distribution with zero mean and a variance given by H13|l . The asymptotic limit of the scaling functions 
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(f>n{mr) will determine the multiscaling of S'2„(r), which are not invariant underi the transformation TIL Hence in 
the OPE of O^p — [?/'(x + r) — V'(x)]^" operators violating TII may appear ^1^. By virtue of TII operators (6^)^" 
(most dominant in the OPE of O^p) do not renormalize and as a result the anomalous dimension of is zero [Toj . 

Therefore, S'^„('') ^ r"^'*"'^' exhibiting only simple scaling as observed in Ref.l?! for 2c?. This conclusion remains 
unaffected even when u is to be given by the solutions of the Navier-Stokes equation. Therefore our results provide 
a symmetry motivated understanding for those in Ref. 6]. It should be noted that in the problem of true passive 
advection, e.g., the advection of a temperature field by an incompressible field, Ofo = and as a result, the condition 
imposed by TII does not exist any more. Therefore, the operators renormalise and contribute to the multiscaling 
of S2n{r) as has been found in Refs.jTH. 

In summary, we have discussed the interplay between the symmetries of the model equations of a binary mixtures and 
its scaling and multiscaling properties in the NESS. We have highlighted the dependence of the statistical properties of 
the NESS on the degree of crosscorrelations between the velocity and gradient fields. Our results open up the intriguing 
possibility of continuously varying multiscaling universality class for such a system. Our results on multiscaling suggest 
the singular nature of the weak scalar limit of the full nonlinear problem and connects the presence of a symmetry 
to such properties. Our results display the key aspects of multiscaling as observed in fluid turbulence. Although the 
noise crosscorrelations (|5)l implies that the velocity field u and the concentration gradient field b have symmetric and 
antisymmetric crosscorrelations u and the concentration field ip have no crosscorrelations. Crosscorrelations between 
u and V' can arise only if fields u are compressible. It should be noted that while comparing with real data one needs 
to impose the irrotational condition with respect to the second index in the third of the Eqs. This arises due 

the fact the vector b is purely solenoidal (b — Wip). This yields a relation between Ng and iV^; however, one of 
them still remains free and allows for variable scaling properties in the steady state. Such properties are essentially 
nonequilibrium since in equilibrium noise crosscorrelations vanish identically. Our results may be verified in turbulent 
experiments with symmetric binary rnixtures in their miscible phase. In particular we suggest experiments active-grid 
generated turbulence experiments 17] with symmetric binary fluid in its miscible phase. Such active-grid experimental 
set-ups, which, to a good extent, are practical realizations of stochastically driven turbulence models, are fitted with 
grids with agitator wings which rotate with frequencies and for durations which are randomly chosen with given 
distributions. Note that in the present case, crosscorrelation functions of the velocity and the concentration gradient 
fields are odd under parity reversal which are absent in equilibrium and can be generated only by appropriate external 
drives. In experimental set-ups such spatial randomness can be incorporated by adding several layers of active grids, 
each layer having random distributions of agitator wings. The effective stochastic forcings, generated by allowing 
agitator wings belonging to different levels of grids to have different rotation speeds and time durations with the 
spatial distributions of agitator wings breaking inversion symmetry will generate non-zero crosscorrelation functions 
of the velocity and the concentration gradient field at the lowest order. In the resulting nonequilibrium steady state, 
the concentration fiuctuations can be studied with light-scattering techniques 1 181 . In numerical experiments such 
crosscorrelations can be generated by extending the methods of A. Sain et al jl9| by controlling the appropriate 
Grashof numbers [23, constructed out of the various noise correlators. In the closing we would like to comment that 
since the ao term in Eq. ^ is 0(1) in symmetric binary mixtures, experimental studies will reveal the active nature 
of the dynamics, and more specifically, measurements of multiscaling should verify that the structure functions >S'2„(r) 
exhibit multiscaling, where as <S'^„(7') exhibit simple scaling even in the weak concentration field limit. Furthermore, 
in asymmetric binary mixtures the correlation length C is small and hence, in the relevant scaling regime over length 
scales ^ C the ao term in Eq. is very small leading to experimental situations revealing passive behaviour. From 
a wider point of view our results illustrate the richness of nonequilibrium phenomenologies. Recently, analyses on 
model nonequlibrium systems revealed the existence of continuous variation of the scaling exponents. Ours is 
an illustration of such variations of multiscaling universality classes in natural systems. Our results suggest further 
studies for similar effects in driven soft-matter and biologically relevant systems. 

The author wishes to thank J. Santos for discussions, D. R. Nelson for critical comments and Alexander von 
Humboldt Foundation, Germany for partial financial support. 
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